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A fluid in the NVT ensemble at T less than the critical temperature T c and 
p = N/V somewhat in excess of Pcoex (density of the saturated gas in the 
gas-liquid transition) is considered. For V — ► oo, a macroscopic liquid droplet 
coexists with surrounding saturated gas according to the lever rule. For fi- 
nite V, droplets can only exist if they exceed a minimum size. A (rounded) 
first order transition of the system occurs when the droplet evaporates into 
the supersaturated gas. Simulation evidence for this transition is given for a 
Lennard- Jones model and interpreted by a phcnomcnological theory. At the 
transition, the chemical potential difference fit — Mcoex scales like 
for a cubic volume V — L d in d dimensions, as L — > oo. 



1.1 Introduction 

There is much interest in the pathway of first order phase transitions, where 
a droplet of the new phase is formed on the background of the old phase 
[1-4] by applying computer simulations [4-9]. For the liquid-gas transition, 
simulating N particles interacting with a Lennard-Jones (LJ) potential in a 
finite volume V and applying periodic boundary conditions in all directions, 
one finds a liquid droplet coexisting with surrounding supersaturated gas for 
an appropriate range of T and p. From such data the excess free energy of 
the droplet and the associated nucleation barrier can be inferred. 

However, for all simulations [10, 11] one must be aware of systematic 
errors due to finite size effects. Nucleation theory considers the (hypothetical) 
unstable equilibrium of one liquid droplet coexisting with an infinite amount 
of surrounding supersaturated gas. This saddle point in configuration space is 
not accessible to a straight-forward simulation, though. For a finite amount of 
surrounding supersaturated gas, the droplet is stable in a suitable parameter 
range [6, 12, 13]. If we increase the total volume (at fixed droplet volume), we 
encounter a size where the droplet is no longer stable but rather evaporating 
(with an increase in supersaturation of the surrounding gas) [6, 12, 13]. 

Some aspects of this transition have recently been proven rigorously [13]. 
Implications are manyfold. The transition is e.g. crucial for the efficiency of 
multicanonical sampling [14] in large systems. It can also help to understand 



2 P. Virnau, L.G. MacDowell, M.Muller, and K. Binder 



related phenomena such as the melting of islands adsorbed on surfaces in the 
submonolayer range [15,16] or microcanonical simulations close to two-phase 
coexistence [17]. 



1.2 Simulation evidence for the Lennard-Jones fluid 

Grand canonical Monte Carlo simulations were carried out in d = 3 dimen- 
sions using a truncated and shifted LJ potential: 
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if r < r c = 2(2 1/6 er) and else. (1.1) 



A multicanonical weighting scheme [18] was applied to simulate states below 
T c . Weighting functions were generated by adopting Wang-Landau sampling 
[19] to our off- lattice ([iVT) ensemble [20]. From the resultant probability 
distribution P^vt(N) one can obtain the chemical potential 



Unvt = -kT[ ^ I + const 



(1.2) 



and A\x = [Invt — Mcoex with /icoex = being an estimate for [i at the 
coexisting curve [21]. Inside the two-phase region and for finite V, A[i jumps 
to a lower value at an effective "spinodal point", which moves towards the 
coexistence curve as V increases [21] (Fig. 1.1). 
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Fig. 1.1. Plot of f3Afi vs. JV, for a LJ fluid in a box of L = 22.5 at T = 0.68 
(An = fiNVT — A'coex). Note that lengths are measured in units of a, temperature 
in units of e/fcs, and /3 = 1/fcsT. o denotes densities at which configurations were 
stored for analysis (compare with Fig. 1.2) 
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Fig. 1.2. Distribution p(N c ) of the cluster size N c for several choices of N (left part) 
and the corresponding distribution of the chemical potential of the supersaturated 
gas p(n) (right part). System parameters are the same as in Fig. 1.1. 



While in principle this has been known for a long time [22], we now present 
direct evidence from simulation. During the multicanonical run, we saved con- 
figurations for typical densities (compare with Fig. 1.1). After simulation, the 
cluster size distribution p(N c ) and the distribution p(p) of the chemical po- 
tential fj, in the gas region outside of the cluster was determined. (Fig. 1.2). 
Any ensemble of atoms whose distance is smaller than 1.5a is assumed to 
belong to the same cluster [23] . The chemical potential was measured by the 
Widom particle insertion method [24]. For N < 350, i.e. on the ascending 
branch of the A\i vs. N curve in Fig. 1.1, p(N c ) is monotonically decreas- 
ing with N c . For N w 355 a peak near N c w 120 appears. It becomes more 
pronounced and moves to larger N c as N increases. This peak represents a 
single large liquid droplet, present in the descending part of the A\x vs. N 
curve. However, the liquid droplet cannot be found in all sampled configu- 
rations: when we sample the distribution of the largest cluster N™ ax in the 
system [21], we find a bimodal distribution with one peak near A^™ ax w 20 
in the whole range of N studied in Fig. 1.2a. This peak corresponds to config- 
urations which consist of supersaturated gas with small clusters but no single 
large droplet. The second peak is identical to the peak of the single droplet as 
shown in Fig. 1.2a. This clear cut evidence for the "evaporation" of the liquid 
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droplet also appears in double peak - distributions of the chemical potential 
(Fig. 1.2b) and the internal energy [21]. 

1.3 Phemonenological theory using the Ising lattice gas 
model and concluding remarks 

The discussion is simplest for the lattice gas, due to its perfect symmetry 
between liquid and gas, brought out in the magnetic interpretation: the two 
coexisting phases have spontaneous magnetization imcoexj and A/i corre- 
sponds to the magnetic field. The thermodynamic potential g s .g(n) of the 
supersaturated gas is simply expanded at the coexistence curve, for magne- 
tization m near -mcoex: 

SS.g(m) = go(T) + ~z ~ (m + mcoex) 2 , (1.3) 

^Xcoex 

Xcoex being the susceptibility (per lattice site) at the coexistence curve. 
Similarly, the state containing a droplet is described by 

SdropM - 9o(T) + ^—(m' + m CO ex) 2 + ^f*/^ , (1-4) 
1 ^Xcoex Li 

where the last term describes the surface free energy of a droplet (per lat- 
tice site) in the capillarity approximation, Sd being the surface area of 
a (i-dimcnsional unit sphere, 7 the interfacial tension, and R* the criti- 
cal droplet radius (satisfying R* = (d — l)Sd"t/[2dmcoex VdH] Vd being 
the volume of the unit sphere), and m' + mcoex = XcocxH for H — > 0. 
The transition occurs for gs.g(m) — 3drop( m ) i^] - R- esca l m g variables 
m = m/ mcoex, H = Hxcoex/ mcoex one finds [13] a first order transition at 
fht = -1 + (d+ l)(^/2d) d /( d + 1 )y7 ( ^ 1)/(d+1) c d /( d + 1 )i- d /( d + 1 \ where the 
field exhibits a jump between — 1 + m t and h[ 2 ^ = (d— 1)H^ /(d+ 1), 
L = L/^ C oex, £coex being the correlation length, and c is a constant that 
is universal near T c {c — 7Xcoex/^coexCcocx [12]}- The variation with 
wag g rs ^. n0 ^ ec j m J5] anc j recently proven rigorously [13]. Of course, 
for finite L this anomalous first order transition is rounded, but it becomes 
sharp as L — > 00 [13]. 
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